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Abstract 

We present a mathematically rigorous quantum-mechanical treatment 
of a one-dimensional nonrelativistic motion of a particle in the potential 
field 

V{x) = gix^^ + g2X^^ , X e R+ = [0, oo) . 

For 32 > and g\ < 0, the potential is known as the Kratzer potential 
Vk{x) and is usually used to describe molecular energy and structure, 
interactions between different molecules, and interactions between non- 
bonded atoms. 

We construct all self-adjoint Schrodinger operators with the poten- 
tial V{x) and represent rigorous solutions of the corresponding spectral 
problems. Solving the first part of the problem, we use a method of spec- 
ifying s.a. extensions by (asymptotic) s.a. boundary conditions. Solving 
spectral problems, we follow the Krein's method of guiding functionals. 
This work is a continuation of our previous works devoted to Coulomb, 
Calogero, and Aharonov-Bohm potentials. 



1 Introduction 

In this article, we present a mathematically rigorous quantum-mechanical (QM) 
treatment of a one-dimensional nonrelativistic motion on a semiaxis of a spinless 
particle of mass m in the potential field 

V{x) = gix-^ + .92a;"^ x eR+ = [0, oo) . (1) 

On the physical level of rigor, the Schrodinger equation with potential ([T|) was 
studied for a long time in connection with different physical problems, see for ex- 
ample [HI [7] and books (TUIIH]- In particular, this potential enters the stationary 
radial Schrodinger equation 



^Pnl (r) = , (2) 
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where n and / are radial and angular quantum numbers, after separating spher- 
ical variables in three-dimensional spherically symmetric QM problems, see e.g. 
[8]. The potential ([ij is singular at the origin, it is repulsive at this point for 
g2 > 0, and has a minimum at a point xq > for 32 > and gi < 0. The 
potential with gi , 52 in the latter range is known as the Kratzer potential [1] . 
The Kratzer potential is conventionally used to describe molecular energy and 
structure, interactions between different molecules [S], and interactions between 
nonbonded atoms [2]. For 32 < and gi > 0, we have the inverse Kratzer 
potential which is conventionally used to describe tunnel effects, scattering of 
charged particles [11] and decays, in particular, molecule ionization and fluores- 
cence [4] . In addition, valence electrons in a hydrogen-like atom are described in 
terms of such a potential When modeling some physical systems, a constant 
is usually added to the angular momentum term, I {I + 1) — f3 + l{l + 1), in 
order to take some effective potential energy into account. For example, in the 
model of a molecule interaction, f3 can represent the dissociation energy of a 
diatomic molecule [5] or, in the scattering problem, this parameter represents 
attractive {(3 < 0) or repulsive {(3 > 0) interactions between charged particles 

m- 

In Figure 1 we show the shape of the potential under consideration for dif- 
ferent values of the parameters. 

Even though a number of works was devoted to the QM problem with the 
potential (H]), a rigorous mathematical analysis of this problem is lacking in 
the literature. The aim of such an analysis (which is, in fact, the aim of the 
present article) is to construct all self-adjoint (s.a. in what follows) Schrodinger 
operators (Hamiltonians) with the potential ([T]) and present rigorous solutions 
of the corresponding spectral problems. 

When solving the first part of the problem, we use a method for specifying 
s.a. differential operators by (asymptotic) s.a. boundary conditions (the so- 
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called alternative method, see [12] )• When solving spectral problems, we follow 
the Krein's method of guiding functional, see [M] and books [15]. This work 
is a continuation of our previous works [5U [3S] devoted to Coulomb, Calogero, 
and Aharonov-Bohm potentials; using the given references, the reader can be- 
come acquainted with necessary basic notions and constructions, like guiding 
functional and Green function. 

As in the above-mentioned works, we start with a s.a. differential operation 
H on R+, 

H = -dl+ gix-^ + .g2a;"^ (3) 

and examining solutions of the corresponding homogeneous differential equation 
{H-W)i! = 0, or 

i)" - {gix-'^ + g2X-^ - W)i}) = Q,W^ |W^|e*'^, < < 27r , (4) 

which is the Schrodinger equation (with omitted factor 2m/ fi?) with a complex 
energy W , for ImVF = 0, we write W = E m what follows. 

The basic operator if+ in associated with H is defined on the nat- 

ural domaiifl -D^(M+) C 

D*^ (M+) = {V^la::) : '(/'*,V'*area.c. inM+; e L2(M+)} , (5) 

it is the adjoint of the so-called initial symmetric operator H associated with H 
and defined on the dense domain Dfj = V (M_|_), the space of smooth functions 
with a compact support, 

V{M.+ ) = {ijj{x) : e C°°(M+), suppxp C [a,/?] c (0,c5o)} (6) 

it is evident that V (R+) C £»^(M+) and H d H+ . The operator H+is gener- 
ally not self-adjoint and even not symmetric; its quadratic asymmetry form is 
denoted by . All possible Hamiltonians associated with H are defined as 
s.a. restrictions of , which simultaneously are s.a. extensions of the symmet- 
ric H , the restrictions to some subspaces (domains) belonging to D*^ (IR+) and 
specified by some additional (asymptotic) s.a. boundary conditions on func- 
tions belonging to D*^ (IR+) under which the asymmetry form becomes 
trivial (vanishes); these domains are maximum subspaces in D*^ {^+) where 
the operator is symmetric!! (see [H]). Our first aim is to describe all these 
Hamiltonians. The special case of gi = corresponds to the Calogero potential 
and was already considered in [21], we therefore keep 7^ in what follows. 

This paper is organized as follows. In sec. [2] we present and discuss some 
exact solutions of equation ([4]) and their asymptotics. In the following five 
sections, we construct all s.a. extensions of H, and perform the corresponding 

^a.c. means absolutely continuous. 

^Although the notions "s.a. extension of and "s.a. restriction of H " are equivalent; 

it is more customary to speak about s.a. extensions; we use one or another of the equivalent 
notions where appropriate. 
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spectral analysis of the Hamiltonians for different ranges of the parameter 32 • 
In sees. (|3.m3.4|) . we consider the case of 32 7^ 0. The special case of 32 = 
is considered in sec. 13.51 In sec. SI we highlight some remarks and possible 
applications of the obtained results. 



2 Exact solutions and asymptotics 

We first consider the Schrodinger equation Introducing a new variable z 
and new functions (t>±{z) instead of the respective x and V' (x), 

z = \x, A = = 2v/^W1e*('^-'')/^ iP{x) = x^/'^^^'e-'/^(|3±{z) , 

„ ^ / V92 + 1/4 , 92 > -1/ 4 , . 

' v/|<?2|-l/4, 52 < -1/4 ' 

we reduce eq. (U) to the confluent hypergeometric equations for (f>±{z), 

zdl(f)±{z) + {I3± - z)d^(j)±{z) - a±(/>±(z) = , 

a± = l/2±^ + 5i/A, /3± = 1±2a* , (8) 

their solutions are the known confluent hypergeometric functions <I'(a±, z) and 
*(a±,^±;^), see [Hill]. 

Solutions V' [x) of eq. (|3]) are restored from solutions of eqs. © by transfor- 
mation (O. In what follows, we use ui{x\ VF), U2{x\ VF), and vi {x; W) defined 

by 

{x; W) = xi/2+A'e--/2<I>(a+, z) = u, (x; W)\^^^^ , 
U2 {x; W) = xi/2-^e--/2$(a_, z) = ^2 (a:; W^)U^-a = ^1 (2^; ' 

V, {x,W) = A2^xV2+Me-/2vi,(«^,^^;,) = \^^^±^u, + ^^2 . (9) 

r(a_) r(a+) 

The function 1*2 is not defined for /3_ = — n, or /x = (n + l)/2,n € Z+, in 
particular, for /i = 1/2. For such ^, we replace by other solutions of eq. (|4]), 
they are considered in the subsequent sections. 

The coefficients of the Taylor expansion of functions u\[x'^W) j x^^'^'^^ and 
U2{x\ W)/x^^^~'^ with respect to x are polynomials in A. Because these functions 
are even in A, the coefficients are polynomials in W, whence it follows that 
Ui {x; W) and U2 (x; W) are entire functions in W at any point x except x — 
for U2 with fi > 1/2. 

If .92 > —1/4 ifJ- > 0), then ui {x; W) and U2 {x; W) are real-entire functions 
of W. If 52 < -1/4 (/X = iM:), then U2 (x; E) = ui (x; E). 

The pairs ui,U2 with /i 7^ Oand for ImVF 7^ are the fundamental 

systems of solutions of eq. ^ because the respective Wronskians are 

Wr {ui,U2) = -2(1, Wr (ui, ui) = -r(/3+)/r(a+) = ~lj{W) . (10) 
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The well-known asymptotics of the special functions $ and see e.g. [TB], 
entering solutions (|9]) allows simply estimating the asymptotic behavior of the 
solutions at the origin, as a; — >■ 0, and at infinity, as a; — >■ oo. 

As a; — 0, we have 

0(a;5/2-M), _i/4 < 32 < 3/4, ^ q, 
7/2(2;; VK)=Ko^^^''"2as(a:) + <( (0</i< 1/2) 

0(x3/2), 32 < -1/4 ifi^lK) 



(11) 



and, if a+ ^ — n, 7^ — m, n, m G 



f ^)x''^-^{l + 0{x)), .92 > 3/4 (a. > 1) 



ui(a;; VK) = < 



where 



-1/4 < 52 < 3/4,52 ^ (0 < < 1, M 1/2) 

"las(^) + ^g^<'/'+*'^"2a.(2;) + 0(a:3/2), 

. 32 < -1/4 (pL = i^f) 

(12) 



wias(x) = (A^o2:)^/'+^ 

(«;oa;)i/2-^ - |^(«o:^)^/2-^ -1/4 < 32 < 3/4, 52 ^ 0, 
W2as(a;) = \ (0 < A* < 1,A* ^ 1/2) 

(^02^)1/2--, 32 < -1/4 (^ = i>^) 



(13a) 



and Kq is an arbitrary, but fixed, parameter of dimension of inverse length. 
As a; — >■ 00, Im 14^ > 0, we have 



Ml 



^ I(^A"+-^+a:»i/V/2(l + 0(.T-i)) = o(2;«el^l'''-'^(^/2)) ^ 



r(a+) 

ui(a;; M^) = \-'^- x-3^'^e-''^{l + 0{x-^)) = o(^-«e-|W'l''' -MW2)) , 
a = 2-i|Ty|-i/23^sin((^/2) . 

The obtained asymptotics are sufficient to allow definite conclusions about 
the deficiency indices of the initial symmetric operator H as functions of the 
parameters (71,52 and thereby about a possible variety of its s.a. extensions. 
It is evident that for Iml4^ > the function ui(a;; W) exponentially increasing 
at infinity and is not square-integrable. The function vi {x; W) exponentially 
decreasing at infinity is not square-integrable at the origin for 52 > 3/4 (^ > 1), 
whereas for 52 < 3/4, it is (moreover, for 52 < 3/4, any solution of eq. Q is 
square-integrable at the origin). Because for ImVF > 0, the functions ui,ui 
form a fundamental system of eq. this equation with ImVF > has no 
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square-integrable solutions for (72 > 3/4, whereas for 32 < 3/4, there exists one 
square- integrable solution, vi {x; W). This means that the deficiency indices of 
the initial symmetric operator H are equal to zero, m± = 0, for (72 > 3/4 and 
are equal to unity, m± = 1, for 52 < 3/4. 

Correspondingly for (72 > 3/4, there is a unique s.a. extension of H, whereas 
for (72 < 3/4, there exists a one-parameter family of s.a. extensions of H. A 
structure of these extensions, in particular, an appearance of their specifying 
asymptotic boundary conditions, depends crucially on a specific range of values 
of the parameter g2- In what follows, we distinguish five such regions and 
consider them separately. 



3 Self-adjoint extensions and spectral analysis 

3.1 The first range ^2 > 3/4 (/i > 1) 

As was mentioned above, the deficiency indices of the initial symmetric operator 
H with (72111 this range are zero. This implies that for g2 > 3/4, the operator 
is s.a. and Hi = H'^ is a unique s.a. extension of H with the domain 

A spectral analysis of the s.a. operator Hi = H^ begins with an evaluation 
of its Green function G {x, y; W) that is the kernel of the integral representation 
of the solution -0* (x) of the inhomogeneous differential equation 

{H - W) V'* (x) = ?7(x) , 7j{x) e (R+ ) (14) 

with Imiy 7^ under the condition that i/'* G D*^ i-C-, that ■0* is square- 

integrabl^, '(/'*(a;) £ L^(R+) (see |24[|25j ). The general solution of this equation 
without the condition of square integrability can be represented as 

ip^x) = aiui(x; W) + 02^1 (x; W) + I{x; W) , 
^:ix)^aiu[ix;W)+a2v'iix;W)+I'ix;W) , (15) 

where 



/*x poo 

I{x; W) = G(+) {x, y; W) v{y)dy + J G^-^ {x, y; W) i^{y)dy , 

pX nOO 

l'{x;W)^ d^G^+'^ {x,y;W)r^iy)dy+ d.,G^-\x,y;W)r^{y)dy , 

Jo Jx 

{x,y; W) = Lu-\W)vi{x; W)ui{y; W) , 
G(-) (x, y; W) = uj-\W)ui{x; W)vi{y; W) , 

with uj given in (jlOp . Using the Cauchy-Bunyakovskii inequality, it is easy 
to show that I{x;W) is bounded as a; — > 00. The condition -0,(2;) G L^(R_(-) 

•^We note, that can be considered as the space of unique square-integrable solu- 

tions of eq. lfT4)l with ImW and any ri(x) G L^(R+). 
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then implies that oi — 0, because ui{x; W) exponentiahy grows while vi{x; W) 
exponentially decreases at infinity. As x — > 0, we have I{x) ~ 0{x^^^), I'{x) ~ 
0(a;^/^) (up to the logarithmic accuracy at 32 — 3/4), whereas vi{x; W) is not 
square- integrable at the origin. The condition € ^^(1^+) then implies that 

02 = 0. In addition, we see that the asymptotic behavior of functions ipt:{x) 
belonging to D*^ (R+) at the origin, as x — )• 0, is estimated by 

M^)^0{x^/'), V:(x)=0(xi/2) . (16) 

Together with the fact that the functions vanish at infinity (see below) , this 
implies that the asymmetry form Ajj+ is trivial, which confirms that in the first 
range the operator is symmetric and therefore self-adjoint (in contrast to 
the next ranges considered in the subsequent sections). 
It follows that the Green's function is given by 

G(x vW)^l G^^H^,y;W),x>y 
G(x,y,i^)-| a(-H^,y;W),x<y ■ 

The representation Q of the function vi in terms of the functions ui and 
U2 is inconvenient sometimes, because the individual summands do not exist 
for some /i although vi does. For our purposes, another representations are 
convenient. For m — 1 < 2/i < m + 1, m > 2, the function vi{x] W) can be 
represented as 

uj(W) 

vi{x; W) = A^{W)ui{x; W) + -^v(^„,){x; W) , 

I (a_) I (a+) 

U(„)(a;; W) = U2 [x; W) - a^{W)T{l3.)u, [x; W) , 

It is easy to see that all the coefficients am{W) are polynomials in W which are 
real for Im W = {W = E). In view of the relation 

hm r(/3)$(a, 13; x) = ^_i_t_Li$(a + n + 1, n + 2; x) 
/3->-n (n + 1)!I [a] 

(see [inillH]), the functions V(m)(a;; W) and Am{W) exist for m — 1 < 2// < to+ 1 
and for any W. In fact, V(^rn) {x] W) are particular solutions of eq. @ which are 
real-entire in W and have the properties (for to— l<2/i<m+l) 

Wr(ui,U(„)) = -2/1, V(^ra){x;W) = x^'^'^'il + 0{x)), x^Q. 

As a guiding functional, wc take 

H^]W)^ U{x]W)£,{x)dx, Dr{M.+ )C^DH, , (17) 
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ui{x;Eo) / V(^^-){y; Eo)^o{y)dy + V(^^){x; Eo) / ui{y; Eo)^o{y)dy 



where U {x; W) — ui {x; W) and -Dr.(K+) is the space of arbitrary functions with 
a support bounded from the right: if {x) S Dr (IR+) supp (p C [0, /3], /3 < oo; 
the domain D is dense in L^(K+). The functional ((T7)l is a simple 

guiding functional, i.e., it satisfies the properties: 1) for a fixed ^, the functional 
$(C;T^) is an entire function of W; 2) if ^{Co;Eq) = 0, Im£:o = 0, S 
then the inhomogeneous equation {H — -Eo)V' = Co has a solution S U; 3) 
4>(i?C; W^) = W^{^;W). It is easy to verify the properties 1) and 3), and it 
remains to verify that the property 2) also holds. Let 

*(Co; Eo) = / ui{x; EMx)dx = 0, e B, supp^o e [0, 6] . (18) 
Jo 

We consider the function ipi^) defined by 

(19)' 

that evidently satisfying the equation {H — Eo)^j{x) — £,o{x). Using condition 
(fT8|) . we obtain that supp?/' G [0, fe], i.e., V G ^r(K+), and therefore, -0 G i^(c, 6) 
for any c > 0. With taking the asymptotic behavior of functions ui(a;;i?o), 
i;(m)(a:; -Eq), and Co(a^) at the origin into account, a simple evaluation of the 
integrals in representation (IT51) gives: 

r 0(xi/2+P)^ 1 < ^ <3 

ip{x)^l 0{x'^/'^\n5), ^i = i ,x^0, 
i 0(x7/2)^ ^ > 3 

i.e., ^/j G Dh,, and therefore, e B. 

The derivative of the spectral function is given by 

a'{E) ^T:-^lm[Lj-^{E + iQ)A^{E + . (20) 

Because uj~^{W)Am{W) is an analytic function of /i, its value at fj, = m/2 
is a limit as/i — m/2. For /i ^ ?Ti/2, representation (l20l) can be simplified to 



^r(a_)r(/3+) ■ 

For S = p2 > 0, p > 0, A = 2pe^"/^ we find 

, /|r(a+)|\^ (2p)2Me-'^9i/2p 

We see that a'{E) is a nonsingular function for > 0. It follows that the 
spectrum of the s.a. Hamiltonian Hi is continuous for all such values of E. 

For E = — < 0, t > 0, A = 2t, the function n{E) is real for all values 
of E where fl{E) is finite, which implies that lmn{E + iO) can differ from zero 
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only at the discrete points En where 1/ ^l{En) = 0. It is easy to see that the 
latter equation is reduced to the equations «+(£'„) = —n,n G Z_|_, which have 
solutions only if 51 < 0, and the solutions En are then given by 

En = + 2^l + 2n)-2 , Tn = \9i\ (1 + 2^ + 2ny^ . (22) 

We thus obtain that for E < 0, the function (j'{E) is equal to zero if 51 > , 
whereas if < 0, this function is given by 



VF^ V^n^A^p n (2r„y'+i r„r(l + 2A. + n) 

^ r(/3+) y |gi|n! 

The final result of this section is as follows. 

For (72 > 3/4 (/i > 1), the spectrum of a unique s.a. operator (Hamiltonian) 
Hi is simple and given by 



speciJi 



l+, gi > 

l+U{En,}, 31 <0 • 



For gi > 0, its generalized eigenfunctions Ue{x) — a' {E)ui{x] E), E > 0, 
form a complete orthonormalizcd system in L^(R+). For 171 < 0, the generalized 
eigenfunctions Ue{x) ~ \J a' [E)u\{x\ E), E > 0, of the continuous spectrum 
and the eigenfunctions C/„(x) — QnUi{x; En), n € Z+, of the discrete spectrum 
form a complete orthonormalizcd system in L^(R_|_). 

3.2 The second range 3/4 > 92 > -1/4, g2 {1 > jj, > 
0, /i ^ 1/2) 

We note that in this section, we consider the range 3/4>52>— 1/4 excluding 
the point g2 — {fi = 1/2), the reason is that the function U2 we use here is not 
defined for ^ — 1/2. The case g2 — {iJ- — 1/2) is considered separately in the 
last subsection. 

The operator /f+ with 172 in the second range is not s.a., and we must 
construct its s.a. reductions. In accordance with the general procedure of the 
alternative method, see [12] and also [25, [25] for examples, we begin with 
evaluating the quadratic asymmetry form Ah+ in terms of quadratic boundary 
forms, which are determined by the asymptotics of functions 'ipsr{x) belonging 
to the natural domain D*^ (K+) at the origin (the left boundary form) and at 
infinity (the right boundary form). Because the potential vanishes at infinity, 
the right boundary form is trivial (zeroj^, see [E], and the asymmetry form 
A/f+ is reduced to (minus) the left boundary form. To determine an asymptotic 
behavior of functions V'* at the origin, we consider these functions as solutions 
of the inhomogeneous eq. (|T4l) with W = Q. Because in the range under 



^Moreover, we can prove that 1/1, vanishes at infinity together with its derivative, 
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consideration, any solution of the homogeneous eq. ^ is square-integrable at 
the origin, the general solution of eq. (fT4|) with W ^ can be represented as 

tp^ix) = aiui{x;0) + 02^2(0;; 0) 

f [Mx; 0)^2(2/; 0) - U2{x; 0)ui{y; 0)] r]{y)dy . (23) 
Jo 

The asymptotic behavior of the functions ui and U2 in representation (j23p 
as X —> is given by ([TT|) and (jlSap . the asymptotic behavior of the integral 
terms is estimated using the Cauchy-Bunyakovskii inequality, and we find 

ip*{x) = aiwias(x) + a2M2as(a;) + 0{x^/^) , 

i^Ux) = a^u[^,{x) + aW^^^ix) + 0{x^'^) . (24) 

With these asymptotics, we calculate the left boundary form [V'*,V'*](0) = 

lima;_s.o(— '!/'*( a^) "01 (a;)+V'*(a;)'0* {x)) and obtain a representation of the quadratic 
asymmetry form as a quadratic form in the coefficients ai and 02 in (j24p : 

= -2/iA:o(ara2 -0201) • 

The coefficients 01,02 are called the (left) asymptotic boundary (a.b.) coeffi- 
cientfU. The requirement on the a.b. coefficients that Aj:/+ vanish results in the 
relatioiH 

02 sin J/ = oi cos , e § (— 7r/2, 7r/2) , (25) 

between these coefficients. It follows that the quadratic asymmetry form 
becomes trivial on the subspaces of D*^ such that the a.b. coefficients of func- 
tions belonging to D*^ satisfy relation (|25|) with fixed v. These subspaces 

are just the domains of s.a. restrictions of , and relation ()25p . with fixed v, 
defines the asymptotic boundary conditions specifying these s.a. operators. 

We thus obtain that for each 52 in the second range, there exists a family 
of s.a. Hamiltonians -ff2,i/ parametrized by the parameter on a circle with the 
domains Dh^^, that are the subspaces of functions belonging to Z)^(R+) and 
having the following asymptotic behavior at the origin, as a; — 0, 

V'(a;) = Ci}''\x) + 0(x3/2) , ^^'{x) = C?/'^'"(a;) + 0{x^''^) , 

%l)'^{x) = uias(fcoa;) sini/ + U2as(a;, fco) cosi/ . (26) 

The spectral analysis of H2,v is similar to that for Hi in the previous section, 
the difference is that the function vi{x]W) is now square-integrable at the 
origin and we must take asymptotic boundary conditions (|26p into account. To 
evaluate the Green's function for H2,u^ we take the representation (fTSj) with 

^The inertia indices of the quadratic form (l/2j/i/to)A_|_ are 1, 1, which confimis the previos 
assertion in sec. ^ that the deficiency indices of H are m± = 1, see I12| . 
^Here and in what follows we use the notation S (a, b) = [a, 6] , a ~ b. 
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ai = for belonging to Dh2^, boundary conditions (P^ . and asymptotics 

(fTTj) . (fT3a|) then yield 



02 



r(2/i) . r(-2Ai)(A/fco) 



2fi 



sm ^ ^ cos i> 



T{a+) r(«_) 
vi{x; W)ri{x)dx . 

Representing the function vi{x; W) in the form 

U2m{x; W) — ky^^^ui{x] W) smv + ky'^^^U2{x] W) cos v , 
U2m{x; W) — —kQ^^^'^ui{x; W) cos v + kl^^~'^U2{x; W) sini/ , 

UJ2 AW) ^uj{W){X/ko)-^smi, + {X/kor^^ cos ly , 

r(a_) 

W2,.(W^) = w(W^)(A/fco)"^cosi/ - (A/fco)^J!7^sini/ , 

r(a_) 

where w is given in (ITOl) . the functions M2,iy(a;; VF) and U2,u{x; W) are real-entire 
in solutions of eq. ^ and U2,i,{x; W) satisfies boundary condition ([26l) . we 
obtain the Green function 

G(a;,y;W^) = (2Aifco)-^r!(W^)M2,.(x; W^)m2,.(2/; VK) 

_l J_ f U2jAx; W)u2,^{y; W), x> y , . 

2^lko\ U2,^{x;W)u2,^.{y;W), X <y ' ^ ' 

where 

n{w) = tJ2,l(w)uj2.AW) . (28) 

We note that the second summand in ([?7)) is real for real W = E. 

As a guiding functional we take the functional $(^; W^) given by (fT7|) with 
/7 (x; 1¥) = U2,,.(a;; M^) and ^ G B = £'r(R+) n fi/a.,- The domain D is dense 
in iy^(R+), D — L^(R+). Following the procedure of the previous section, we 
show that $(^; z) is a simple guiding functional, i.e., satisfies the properties 1)- 
3) cited in subsec. 13.11 It is easy to verify the properties 1) and 3). We prove 
that the property 2) also holds. Let 

^o-Eo) - / U2Ax;Eo)UAdx = 0, ^0 e D, supp^o e [0,b] . (29) 
We consider the function 



tjj{x'' 



2/iKo 



U2,u{x;Eo) I U2,u{y; Eo)^o{y)dy + U2,u{x; Ea) I U2,,y{y; Eo)^o{y)dy 
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which is a sohition of equation 

[H - E^)i;(x) ^ . 

Using condition (I^H]) . we obtain that supp^ G [0, &], i.e., ij] E -Dr(K+), and 
therefore G L^{c, b) for any c > 0. 

The function ip{x) allows the representation 



-U2,^(x;Eo) [ U2,„{y;Eo)^o{y)dy, c= [ U2,u{y; Eo)^n{y)dy . (30) 

Jo Z^Kq Jq 

Using the asymptotics of functions U2,,y{x; Eq), U2,i/(a;; Eq), and ^0(2;) and simple 
estimates of the asymptotic behavior of the integral terms at the origin, we 
obtain that the asymptotic of ipi^) at the origin is given by 

i^{x) = cw2,.(x; En) + ©(x^/^-''), a; ^ , 

which implies that ip S Dh^ „ and therefore -0 g D. 
The derivative of the spectral function reads 

a'{E) = {27rfikoy^lmn{E + iO) . 

It is convenient to consider the cases < tt/2 and v — ±7r/2 separately. 
We first consider the case v = tt/2 where we have 

u2.^/2{x;W) kl^^^^ui{x;W), 

For = p2 > 0, p > 0, A = 2pe^*''/^ we have 

|r(a+)|2 (2p/fco)2^e-'^fi/2p 



a'{E) 



r2(/3+) 27rfco 



such that (j'{E) is finite and specH^^/^ = M+. 

For E = — < 0, t > 0, A = 2t, the function is real for all values of E 

where U{E) is finite, which implies that lmn{E + iO) can differ from zero only 
at the discrete points £"„ where 1/ fl{En) — 0. The latter equation is reduced 
to the equations «_(£"„) = —n,n G (r(a-|_) = 00) or 

l + 2n + gi/Tn = -2n, neZ+ . (31) 

Eqs. ([3T|) have no solutions for > and for < we have (we will denote 
the points of discrete spectrum for ly — ±tt/2 by £"„) 

I51I r. _2_ gf 

— ~n 



l + 2^j + 2n' " (l + 2/i + 2n)2 



12 



such that we obtain 



^ (^) - ^^^"^^^ - ' ^" fC^:^ V (l + 2M + 2n)n! " 

It is easy to see that for the case of i/ = — 7r/2, we obtain the same results 
for spectrum and eigenfunctions as it must be. 

The final result for the Hamiltonian if2,±7r/2 is as follows. Its spectrum is 
simple and given by 

fj. _ ( R+, gi > , 
specii2,±./2 - I u {£„, n e Z+}, g, < ' 

and the complete orthonormalized system of its eigenfunctions in L^(IR+) is 
given by 

Ue (x) = ^a'{E)kl'^+^ui{x; E), E>0 , 

2l-M|£ \3/i-l^/2 

_r r-i(i + ")r(i + "-2/i), o<M<i/2 

^" " \ r-i(2 + n)r(2 + 71- 2m), 1/2 < /X < 1 ' ^ ^+ ' 
for gi > 0, and by 

Ue {x) = ^/^^ky^+''ui{x;E), E>0, 

for gi < 0. 

Now, we turn to the case < 7r/2. In this case we have 
a'{E) = {2n i^ko cos^ ImF^ I {E + iO) , 
F.AW) = HW) + tan. UW) = '•'^-'^^-''^f • 

r(^)r(a_) 

For = p2 > 0, p > 0, A = 2pe-^''/'^, we have 

^'^^ " 27rfco cos2 jy[A2(£;) + ^^b^^e)] ' ^^^^ 

where = ReF2,v{E) and AtS(£') = - Imi^2,i/(-E'). A direct calculation 

gives 

= ^^^StS^ (e-".-^cos(.,., + e™/.) , 
^^^^ jrK)PWe-.W.. ^^ (33) 
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For E = — < 0, r > 0, A = 2t, the function F2^^{E) is real, there- 
fore, cr'{E) can differ from zero only at the discrete points En{v) such that 
F2,v{En{v)) = 0, or f2{En{v)) = — tani/, and we obtain that (derivatives with 
respect to E are denoted by primes in eq. 

a'{E) - [-^t^koF^jEnM) cos^ SiE - E^{v)) , 

n 

HAEnH) = f2iEnil^)) < 0, dMl,) = - COS-2 [/s (^^n (^) )] " ' > • (34) 

I. Let gi > 

For = p2 > 0, p > 0, the function cr'(£') ([321) is a finite positive function. 
At £; = 0, wc have B{0) = and 

mUuo = 0, tanj.0 = -r(/3_)(5i/fco)'^r-i(/3+) . 
It is easy to see that 

F2AW) = tani/- tani^o - (2^fco cos^ i^o) + O(W^), , 

3i(5i/fco)-^ / 3r(l + 2At) 



/icosi/Q y 2fco(i + 2/i)r(2-2/x) ■ 

It follows that for ly i'q, the function cr'{E) is finite at £' = 0. But for — i^o 
and for small -E, we have: 

a'{E) = -i^-^ Im (£; + iOy^ + 0(1) = 'f^SiE) + 0(1) , 

TT 

which means that there is the eigenvalue = in the spectrum of the s.a. 
Hamiltonian H2,uo ■ 

For E = -t2 < 0, a = 2t, the function f2{E), 

f (F\ = r(/3,)r(l/2 + M + gi/2r)(2r/fco)^^ 
^ r(/?+) r(l/2-/i + 5i/2r) 

has the properties: f2{E) is smooth function for E e (— oo,0), f2{E) ^ 00 as 
E = —00, /2(0) = — tani^o- Because f2{E^\n) < 0, see eq. (p4l) . the straight 
line /(-E) = 2^iai\v,E G (— oo,0], can intersect the plot of the function f2{E) 
no more than once. 

That is why the equation F2^,j{E) = has no solutions for v G (i/o,7r/2) 
while for any fixed v £ (— 7r/2, z^oli this equation has only one solution (i/) £ 
(—00, 0], which increases monotonically from —00 to as changes from —t:/2 + 

to l-Q. 

We thus obtain that the spectrum of H2,iy, W\ < 7r/2, with gi > is simple 
and given by 

SDCci?, - / I»+ U {E^-^ M}, ^ e (-^/2, 

specH2,, - I ^ ^ (i^o, 7r/2) or ly = ±n/2 ' ^^^^ 
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The generalized eigenfunctions 



UEix) = ^a'{E)u2Ax]E), E>0 , 
and (for u e (— 7r/2, vq]) the eigenfunction 

Unix) = U{x) = \^-2tikoF^jE'^-^ (i^))) cos^ z.] U2,. (a;; 

of i?2,i/, form a complete orthonormalized systems in L^(R+). 
II. Let 51 < 0. Then: 

For E ^ > 0, p > 0, X ^ 2pe-*^/^ formulas ([321) and §^ hold true. 
Because the functions A{E) and are finite at = (-B(O) 7^ 0), the 

function cr' [E) ([5^ is a finite positive function for E >Q. This means that for 
E the spectra of s.a. Hamiltonians H2.V are simple, purely continuous, and 
given by speci?2,iy = 

For E = -T^ < 0, T > 0, A = 2t, we have 

f ^ r(/3_)r(l/2 + A^-|.gi|/2r)(2T/fco)2A. 
^ r(/3+) r(l/2-At-|.gi|/2r) 

It is easy to see that for fixed v, the spectrum is bounded from below and 
the equation F2^u{E) = has infinite number of solutions 

En [v) = -g\l^r? + 0(n-3) , (36) 



asymptotically coinciding with p2|) as n — 00. 

We thus obtain that the spectrum of H2,v, \v\ < 7r/2, with < is sim- 
ple and given by specH2,v = U {En (j^)}- The corresponding generalized 
eigenfunctions of the continuous spectrum 



Ueix) - yJa'{E)u2Ax]E), E>0 , 
and eigenfunctions of the discrete spectrum 

Unix) = [-2fikoF^jEn (ly)) cos^ U2,.(x; En (l^)). En H < , 

of /f2,i/ form a complete orthonormalized system in L^(]R+). 

It is possible to give a comparison description of the Hamiltonians H2.1,, < 
7r/2 in more detail. 

The function f2{E) has the properties: f2{E) 00 a.s E ~oo; /2 (£„ ± 0) = 
±(xi, n G Taking the third equality in p4p into account, we can see that: 
in each energy interval {£n-i, £n), n G for a fixed S (— 7r/2, 7r/2), there 
are one discrete level En{v) which increases monotonically from £n-i + to 
f „ — when v changes from 7r/2 — to — 7r/2 + (we set £_i = —00). We note 
that the relations 

hm En [v) = lim En+i {v) = £n, n e Z+ , 

I/— >7r/2 1^— > — 7r/2 
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confirm the equivalence of s.a. extensions with parameters v = —tt/2 and 
V = 7r/2. 

We note that it is possible to find the explicit expressions for spectrum, 
spectral function, and the complete orthonormalized system of (generalized) 
functions of the s.a. Hamiltonian for i/ = 0. In this case, results are the same 
as in the first range (52 > 3/4) with additional change fi — /i. One can easily 
verify that such calculated spectrum coincide with the spectrum {En (0)}. 

It should be also pointed out that bound states exist even for the repulsive 
potential, 52,51 > 0, see the dashed line on the Figure 1. 



3.3 The third range g2 = -1/4 (/x = 0) 

The analysis in this section is similar to that in the previous one, a peculiarity 
is that a+ = a_ = a = 1/2 + 51/A, /3+ = /3_ = 1, ui (cc; W) = U2 {x; W), 
and representation ^ of vi {x; W) in terms of ui and 1*2 does not hold. As the 
solutions of eq. (HJ with 52 = —1/4, we therefore use the functions ui{x;W), 
U3{x; W), and vi {x; W) respectively defined by 



U3 [x; W) = x^^^e-'/^— [:e^$(1/2 + fi + .gi/A, 1 + 2fi; z)]^^„ + m {x; W) In ko 



ii {x; W) = x^/^e-'/H{a, 1; z) = m {x; W)\^^_^ , 
d_ 

vi (x; W) = x^^^e-'^/Hia, 1; z) = r-'^{a) \iOa{W)uf'> (x; W) ~ U3 (x; W) 
wo(VF) = 2^(1) - ^(a) - ln(A/fco), a - 1/2 + gi/X , 

where ^p{a) — V (a)/T{a) and fco is a constant. The functions ui {x;W) and 
U3 (x; W) are real entire in W. 

The asymptotic behavior of these functions at the origin and at infinity is 
respectively as follows. 

As X -> 0, z = Ax — > 0, we have 

ui(x; W) - fc(7'/\ias(a:) + 0{x^/^), ui^x) = (kox)^/^ , 

U3(x; W) = fc(7'/'u3as(x) + 0(x3/2 ^ (fcox)l/2 ^ 

vi{x; W) = k^'/^T-'ia) [L^o(M^)wias (x) - wgas (x)] + 0(x3/2 Inx) . (37) 

As X — > 00, Im W > 0, we have 

Mi(x; W) = r'i(a)A""^xSi/^e^/2 [l + O(x-i)] ^ 00 , 

vi{x; W) = A-"x-fi/-^e-^/2 _^ O(x-i)] ^ . (38) 

The functions ui and U3 are linearly independent and form a fundamental 
system of solutions of eq. (jlj, as well as the functions ui and vi for ImVF ^ 0, 
see sec. [2j 

Wr(ui,?/3) = 1, Wr(ui,ui) = -r-i(a) . 
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We recall that, for g2 = —1/4, the deficiency indices of the initial symmetric 
operator H are m± — 1, and therefore there exists a one-parameter family of 
s.a. extensions of H with (72 = —1/4; see sec. [H 

To evaluate the asymmetry form in terms of a.b. coefficients, we need to 
determine the asymptotics of functions tp^ belonging to the natural domain 
£'^(M+) at the origin. To this end, we use representation (1^51) of the general 
solution of eq. (IT4l) with W = Q where the natural substitutions 021/2 — 02^3 
and M2/2// — >■ — M3 must be made. Using the Cauchy-Bunyakovskii inequality for 
estimating the integral terms, we obtain that the desired asymptotic as a; —J' 
is given by 

V'*(a;) = aiMias(a;) + a2U3as(a;) + 0(a;^^^ Inx) , 

and we fincfl A^+('0*) — A:o(aTa2 — 0201), the coefficients 01,02 are just a.b. 
coefficients. The requirement that vanish results in the relation 

flicosi^ = a2Sini?, e § (-7r/2, 7r/2) . 

This relation with fixed •& defines the domain of a possible Hamiltonian as a s.a. 
restriction of or a s.a. extension of H . 

The final result is that: for 32 = — 1/4 , there exists a family of s.a. Hamil- 
tonians i/3,,3 with the domains 

Dh^.^ = {V' : V' e £'^(M+), iP satisfies ^) , 

where l|39p are the asymptotic s.a. boundary conditions at the origin 

^ = C4>3,^^s(.x) + Oix^^'' Inx), ^' = CV'a^^asl^) + 0{x^'^ Inx) a; ^ , 
■03,i?as(a;) = Mias(a;) sin -9 + U3as(a;) cos d . (39) 

To evaluate the Green's function G{x,y; W) for Hs^^, we take the represen- 
tation (fT5|) with ai — for belonging to Dh^^ C -D^(M+), boundary 
conditions (p9| and asymptotics (|37|) then yield 

/•oo 

a2 -r^(Q!)cosi9[a;o(VK)cos?9 + sin?9]"^ / vi{x;W)ri{x)dx . 

Jo 

Using the representation 

T{a)vi — (wq sin — cos i9)u3.ij -I- (wq cos •& + sin i?)'U3_tf , 
"3,1? (x; W) = ui{x;W) sin-d + U3{x;W) cos -d , 
W3,fl(a;; VT) = ui{x; W) cos'd — 1/3(2;; W) sind , 

^This structure of Ajj+ confirms tfie previos assertion tliat the deficiency indices of H are 
m± = 1. 
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where u^^s and u^^^ are solutions of eq. Q real-entire in W, and u^^^ satisfies 
boundary condition ((39)) . we find 



G{x, y- W) = Vl{W)u^^^{x; W)u^,^{y; W) 

u^^i){x;W)u^^^{y;W), x>y 
U3^i){x;W)u3^^{y;W), x<y ' 

n{W) = {ujo{W)cos^ + smd)-\ujQ{W)smd-cosd) . (40) 

We note that the second summand in G{x, y; W) is real for real W — E. 

It is easy to verify that the guiding functional given by p7)) with U = ms,^ 
satisfies the properties 1) and 3) cited in subsec. 13.11 The proof that it satisfies 
the property 2) is identical to that presented in subsec. 13.21 for the second range 
1 > /X > 0. It follows that the spectra of H^^^ are simple. 

The derivative of the spectral function is given by cr'{E) = tt^^ Im [^1{E + iO)]. 

We first consider the case d = tt/2 where we have 

UZ,T,/2{X]W) = ui{x]W) , 

a'{E) = -T:-^lmVL{E + n{W) = V^a) + ln(A/fco) . 
For = p2 > 0^ p > 0, A = 2pe~*^/^ we find 

.'(i.)4(l-tanh!^)>0. 

For E = -r^ < 0, r > 0, A = 2t, and gi > 0, the function n{E) is of the 
form 

n{E) ^ i;{l/2 + gi/2r) + ln(2T/fco) , 

which implies that for gi > 0, there is no negative part of the spectrum. 
For E = -T^ < 0, r > 0, A = 2t, and gi < 0, we have 

n{E) = ^(1/2 - |.gi|/2r) + ln(2r/fco), lmn{E) = Im V(l/2 - |gi|/2r) , 

which implies that there are discrete negative energy levels in the spectrum, 

Sn = -.9?(1 + 2n)-^ r„ = \gi\{l + 2n)-\ n e Z+ 
<j'{E) - Qn5{E ~ £n), Qn = 2\g,\ (1 + 2n)-3/2 _ 

It is easy to see that for the case of = ~tt/2, we obtain the same results 
for spectrum and eigenfunctions as it must be. 

We thus obtain that for gi > 0, the spectrum of i?3,±-7r/2 is simple, continu- 
ous, and given by specif3_±7r/2 = and a complete orthonormalized system 
in I/^(IR+) of its generalized eigenfunctions consists of functions 



Ue{x) = y/a'{E)ui{x;E), E>0 
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For gi < 0, the spectrum of J?3,±7r/2 is simple and given by species. ±7r/2 = 
K+ U {£n, n G and a complete orthonormalized system in L^(M+) of its 

(generalized) eigenfunctions consists of functions 

Ue{x) = ^/a'{E)ui(x;E), E>0 , 

Unix) = 2\gi\ (1 + 2n)"'/' ui(x;£„), £„ < . 

We note that the spectrum and eigenfunctions for coincide with those 

for He with (?2 > 3/4, if we set /i = in the respective formulas in subsec. 13.11 
We now turn to the case < 7t/2. In this case, <j'{E) can be represented 

as 

a'{E) = (7rcos2z?)-ilm[w3(i; + i0)]"^ , 
W3(VK) = V(a) + ln(A/fco) - 2V^(1) - tant? . 

For = p2 > 0, p > 0, A = 2pe-*''/^ and gi < 0, we have 

B(^E) 

" TTCOS^ ^[A^{E)+B^E)] ' ^^^^ 
where uj^lE) — A{E) — iB{E). The function B{E) can be explicitly calculated: 

B{E) = ^ ( 1 - tanh ) > 0, V£; > , (42) 



2 V 2^/E 

whence it follows that for all E >Q, the spectrum of Hj,^^ is purely continuous. 

For = > 0, p > 0, A = 2pe~^'"/'^, and gi > 0, the spectral function is 
given by the same eqs. (|41l) and l^^ . But in this case, B{Q) = and the limit 
\\mv/ ijJz{W) must be carefully examined. 

At small W , we have 

UJ3{W) = (tan?9o - tani?) - (6.gJ)"^ + 0(1^^), taniJo = ln(gi/fco) - 2?/>(l) . 

For ^ z^o, the function (t'{E) is finite at i? = 0. But for d = -Oq and small 
we have 

<^'{E) = ^ Im {E + zO)-^ + 0(1) = -^<5(£;) + 0(1) , 

TT COS^ Wo COS^ Wo 

which means that the spectrum of the Hamiltonian H^^^g contains an eigenvalue 
E = 0. 

For E — — < 0, T > 0, \ — 2t, the function ujsiE) is real, therefore, (t'{E) 
can differ from zero only at zero-points £'„ = £'„ (i9) of W3(i?) (a;3(£'„) = 0), 
which yields 

a' (E) =Y, [-ko^'3iEn)cos^ ^]~^ S{E ~ En), cj^(-B„) < , 

n 

di,En (i?) = [cos^ z?w^(S„)] < . (43) 
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For gi > 0, we have 

uj;j{E) = V(l/2 + .gi/2r)+ln(2r/5i)+tani9o-tanz? , 
UJ3{E) = (l/2)ln|£:| - tan?9 + 0(l), £^^-00, 
^^3(0) = tan')?o — tan-i? . 

For d < -do, the equation uj3{E) = has no solution, whereas for d > -do, it 
has only one solution £;(-) (1?). Because eq. 1^ holds for di,E^-^ {1}), E'--^ (1?) 
increases from —00 to when d changes from 7r/2 — to i^o- 

For (71 < 0, we have 

uJsiE) ^ V(l/2 " |<?i|/2t) + ln(2T/fco) - 2tjj{l) -tan-d , 
uj^{E) = (1/2) In |£;| - tan 7? + 0(1), £: ^ -00 . 

It is easy to verify that the equation co^^E) = has an infinite number of 
solutions En , n e Z+ , bounded from below and asymptotically coinciding with 
^ as n 00, En = -5?/4n2 + 0{n-^). 

We thus obtain that for gi > 0, the spectrum of H^^^ is simple and given by 
speciJa^^ — U {E''^^ and a complete orthonormalized system in 
of its (generalized) eigenfunctions consists of functions 



Ue{x) - ^a'{E)u3^4x;E), E>0 
U(x)-- 



A:oCos^9^^(S(-) i^)) 



(the eigenvalue (i3) exists, and therefore E'^^'^ {■&) and the corresponding 
eigenfunction U{x) enter the inversion formulas only if > i?o); for gi < 0, the 
spectrum of H^^^ is simple and given by spec7?3_^ = K_|_ U and a complete 

orthonormalized system in L^(R+) of its (generalized) eigenfunctions consists 
of functions 

Ue{x) - ^<j'{E)u3Ax;E), E>0 , 

Unix) ^ [-koCOS^duj'^iEn)] Us^XlEn), En < . 

It is possible to describe the discrete spectrum for < 7r/2 and < in 
more details. To this end, we represent the equation Lu^iE) = in the equivalent 
form 

f^iE) = tan 7?, h{E) = ^(1/2 - \g,\/2T) + ln(2r/fco) - 2V'(1) . 
Then we have 

/(-oo) = cx), / (Sn ± 0) = ±oo, n e Z+. 

Because eq. (|43| holds, we can see that in each interval (f„, £n+i), n G {^1} U 
Z+, there is one discrete eigenvalue En and En increases monotonically from 
f„ + to £„+i — when § changes from 7r/2 — to — •7r/2 + (we set f_i = — oo). 
We note the relations 

lim En-ii-d) = lim S„(z9) . 
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3.4 The fourth range §2 < —1/4 (/i = ix, x > 0) 

The analysis in this section is completely similar to that in Section !??^ (although 
the results for the spectrum differ drastically). We therefore briefly outline basic 
points. 

According to Section ^ , the deficiency indices of the initial symmetric oper- 
ator H with g2 < —1/4 are m± = 1, and therefore there exists a one-parameter 
family of its s.a. extensions. 

To evaluate the asymmetry form , we determine the asymptotics of 
functions "0* belonging to Z)^(R+) at the origin using representation (|23|) with 
(1 — i>i oi the general solution of eq. (ITil) with W = Q and estimating the 
integral terms by means of the Cauchy-Bunyakovskii inequality, which yields 

i'^.ix) = aiui^{x) + a2U2as{x) + 0(a;^/^), a; -> , 
V-Ux) = aiu\^^{x) + a2M2as(a^) + 0(a;^/2), a; ^- , 
uiUx) = (fcoa;) W2as(x) = (fcoa:)^/'-'" - ^I^I^R , (44) 

and we fined ii'*) — — 2i>r(aiai —0202). The requirement that A/j+ vanish 
results in the relation oi — e^*^a2, G S (0, tt) defining the domains of possible 
s.a. Hamiltonians. 

The final result is that: for each 52 in the range 52 < ^1/4, there exists a 
family of s.a. Hamiltonians H^^g with the domains 

Dh^o = {V' : e D*^{R+), V satisfies (gS])} , 

where (j45p are the asymptotic s.a. boundary conditions at the origin 

= CV'4as(a;) + 0(x3/2), ^ c<,3(x) + 0(x1/2), X ^ , 



'4'4as{x) = e* Ulas(a;) -|- C ' M2as(a;) = ?/'4as(x) . (45) 

To evaluate the Green's function G{x,y;W) for Hi^e, we use representation 
with fli = for tp^,(x) belonging to Dj^^ g C 15*^(1^+), boundary conditions 
(|26|) and asymptotics (|44p then yield 



^ °^ ' vi{x;W)7]{x)dx, Ui^eiW) ^a{W) + b{W) 



uj{W)loa,0{W) Jo 
am = e-mWMr:, ,(^) _ „_.r(/3_)(A/fco)- 



r(a) ' ' ' r(a_) 

Using the representation 

vi{x; W) = _ (Vfco)'"fco ^ [iw4,e(VK)u4,e(a;; VK) + w4,e(M^)M4,e(x; W^)] , 

4>f 



^This structure of Ajj+ confirms that the deficiency indices of H are m± = 1. 
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where 



<I;4 g(t^) = a{W) - b{W), U4,e(^; W) = i\(r''^k]l^-'''u^{x- W) - e''^ k^^^+^'^uiix; 

where U4_e and U4,e are solutions of eq. ^ real-entire in W , and u/^fi satisfies 
boundary conditions ()45p . we find 



4>tffco \ U4,e(a;; W^)w4,e(2/; W^), x<y ' AxkoWi^eiW) ' 

the second summand in G{x, y; W) is real for real W = E. 

It is easy to verify that the guiding functional given by (1171) with U = U4^g 
satisfies the properties 1)- 3) cited in subsec. 13. 1[ whence it follows that the 
spectra of Hi^e are simple. 

The derivative of the spectral function is given by a'{E) = tt^^ IniJ7(£' + iO). 

For E = p'^ >0,p>0, X = 2pe-*''/2^ ^nd gi < 0, we have 

, ^ {4Tr>ikoyUl-\D{E)\^) , ^ 

a'(E) = TT-^ ImniE) = ' _ ^' ^ , (46) 

^ ' ^ ' (1 + D{E)){1 + D{E)) ^ ^ 

g-2»0p(^)p(^_)g2i^ln(feo/2p)g-,r>. 

D{E) = - 



biE) r(/3_)r(a) 
Because 

1 _L p-27rx-p-7rgi/p 

\D{E)\^^—- _<1, Vp>0, (47) 

specH^ g = R+ and is simple. 

For = p2 > 0, p > 0, A = 2pe"*''/2, and gi > expressions (gSl) and (jiT)) 
for <j'{E) hold true. But in this case, we have |£'(0)| = 1 and must carefully 
examine the limit \i'a\\Y^Q^l{W). 

It is easy to see that for small W , we have the representation 

■ 1 I e2*(eo-e) o 2 



4>f/co [1 - e2i(»o-e)] + iW/A ' " + Ax"^) 



2z r(/3_) ' 

where [(/j/tt] is the entire part of (p/n. For 6 7^ 6*0, the function <j'{E) is finite at 
E^Q. But for 61 = Oq, we find 

cr'(£; + 0) = -TT-i (A/2><fco) Im (£; + iQY^ + 0(1) = {A/2xka) 5{E) + 0(1) , 

which means that the spectrum of the Hamiltonian H/^^^ with gi > contains 
the eigenvalue E = Q. 
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For E = —T^ < 0, T > 0, A = 2t, the function fl can be represented as 
n{E) = 7rtane(£), e(£;) = 6 + er- er{E) + xln(fco/2T) , 

where 

er{E) = ^[lnr(l/2 + .gi/2T + i><) -lnr(l/2 + .gi/2r-i>f)] 

r -7r|5i|/2T + xln(|5i|/2T) + 0(l), 51 <0 
\ xln(3i/2T) + 0(T), 51 >0 , ii^u ^ 

^r(-oo) = ^ In nig±g + 0(1/t), £ ^ -oo 
The asymptotic behavior of 9(£') at the origin and at minus infinity is given by 
" w\gi\/2T + 0{l), 51 <0 



e{E) = { \ e + er + K\n{ko/9i) + 0{T), gi>0 ' . 

9 + Or- evi-oo) + >fln(A:o/2r) + 0(l/r), £; ^ -oo 

Because U,{E) is a real function for E < a' (E) can differ from zero only 
at the points En = En{9) where 0(i?„) = n/2 + im, n G Z, which yields 

cj\E) =Y,QlKE - En), Qn = [4xfcoe'(£;„)]"'/' , Q'{En) > 0. 

n 

We can obtain an additional information about the discrete spectrum of 
H^^g. Representing the equation 6(£'„) = ■jt/2 + irn, n G Z, in the equivalent 
form 

fiiEn) = Tr/2 + n{n - 9/tt), U{E) = 0r - OriE) + xln{ko/2T) , 
deEniO) = - [/i(i?„W)]"' = - [Q\En{e))]-^ < , 
we can see that the following assertions hold. 

a) The eigenvalue En{0) with fixed n decreases monotonically from En (0) 
to En (tt) — when 9 changes from to tt — 0. In particular, we have 
En-i{e)<En{e),'in. 

b) For any 51, the spectrum is unbounded from below: En — > —00 as n — >■ 
—00. 

c) For any 9, the negative part of the spectrum is of the form En = — fcgm^e^'^l"!/* 
0(l/n)) as n —00, where m = m(5i, 52, 9) is a scale factor, and asymp- 
totically (as n —00) coincides with the negative part of the spectrum 

in the Calogero model with coupling constant 52 under an appropriate 
identification of scale factors. 

d) For gi < 0, the discrete part of the spectrum has an accumulation point 
£■ = 0. More specifically, the spectrum is of the form En = —gf/An^ + 
0{l/n^) as n ^ 00 (as in all the previous ranges of the parameter 52) and 
asymptotically coincides with the spectrum for 52 = 0, see below. 
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e) For gi > 0, the discrete spectrum has no finite accumulation points. In 
particular, possible values of n are restricted from above, n < rimax, where 

_ / /4(0)/7r - 1/2 if /4(0)/7r - 1/2 is integer 

- [ [/4(0)/7r + 1/2] if /4(0)/^ - 1/2 > [/4(0)/7r - 1/2] ' 

and the level E = is present in the spectrum for 6 ^ 6q only. 

The final result is as follows: the spectrum of H4,g is simple and given by 
specH4^e = R+ U {£"„ < 0}, — oo < n < rima,^, where rimax < oo for 51 > 
and Umax = 00 for (71 < 0, and the set of the corresponding (generalized) 
eigenfunctions 

Ue{x) - y/a'{E)ui^g{x;E), E > 0; Unix) = Q„'a4,e(x; E^ < 0, 
form a complete orthonormalized system in 

3.5 The fifth range ^2 = (/i = 1/2) 

The analysis in this section is similar to that in subsec. 13.21 A peculiarity is that 
the function U2 is not defined for fi = 1/2, and we therefore use the following 
solutions of eq. (H)): 

ui{x;W) = xe~^/'^<i>{ai/2,2;z), U5{x;W) = ■a5{x;W) - gi\nkoui(x;W) , 
viix; W) = xe~^/2*(ai/2, 2; z) = T-^a,/^) [0.1/2(^^)^1(2;; W) + ur,{x; W)] , 

where 

ai/2 = 1 + gi/X, 

U5{x; W) = e-^/2xi/2 [x-^$(a_, z) + g,T{p^)x'^^a+, z)] ^^^^^ , 
^1/2(1^) = .giC+gi [ipiai/2) + ln(A/fco)] - gi - A/2 , 

C is the Euler constant. The asymptotics of these functions at the origin and 
at infinity are respectively as follows. 
As X — >■ 0, z = Ax — >■ 0, we have 

ui{x; W) = kg'^uii^{x) + 0{x'^), U5{x;W) = U5as(2:) + 0{x'^ lux) , 
vi{x; W) = r-i(ai/2) [k^^oji/2{W)ui^s{x) + U5as(a;)] + 0{x^ Inx) , 
uinsix) = kox, M5as(a;) = 1 + gixln{kox) + Cgix . (48) 

As a; — ?► oo, Im ly > 0, we have 

mix; W) = r-i(ai/2)A-i+5i/^a;+^'i/-^e^/2(l + 0{x-^)) oo , 
viix; W) = A-^'i/-^a;-fi/^e-^/2(^ _^ 0{x-^)) ^ . 
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The functions ui {x; W) and U5 (x; W) are real-entire in W. These functions 
form a fundamental system of solutions of eq. the same holds for the 

functions ui,vi for ImVF ^ 0, see subsec. 



Wr(ui,U5) = -1, Wr(Mi,Ui) = -l/r(ai/2) = -w(VK) . 

As we know from subsec. 13.21 for (72 < —1/4, the deficiency indices of the 
initial symmetric operator H are m± = 1, and therefore there exists a one- 
parameter family of its s.a. extensions. 

For evaluating the asymmetry form , we determine the asymptotics of 
functions tp^,, belonging to D*^{R+), at the origin using representation (^51) of 
the general solution of eq. (fn]) with W = 0, where the natural substitutions 
02^2 — > a2ii5 and U2/2fi — must be made, and estimating the integral terms 
by means of the Cauchy-Bunyakovskii inequality, which yields 

ip^ix) = ai-uias(a;) + a2U5as(a;) + 0(x^/^) , 

= ai<,(x) + a2u'^,,{x) + 0{x^^^) , (49) 

and we fincj^ — — fco(aTa2 — 0201). The requirement that vanish 

results in the relation oicose = 02 sine, e S S (— 7r/2, 7r/2) . 

The final result is that for 172 = 0, there exists a family of s.a. Hamiltonians 
H^^t with the domains 

Dh,_, = {ij-.-ipG D*^{^+), ^ satisfies ^) , 
where (jSOp are the asymptotic s.a. boundary conditions at the origin 

^ - CV5,eas(a:) + 0(x3/2), ^, ^ c^'^^^,,{x) + 0{x^'^), x^Q, 

i>5.easix) = Uiiisikox) Sme + U5as{x) COS £ . (50) 



To find the Green's function G{x,y; W) for i/s^e, we use representation 
with ai = for ^^{x) belonging to Dn^^f, C D*^(M+), boundary conditions ([50 
and asymptotics (|48| then yield 



r^(ai/2)cose 

02 = -— ; : — / viix-^Wj-qixJdx 

i^i/2(yy)cose - ko sine Jq 



Using the representation 



ko^{ai/2)vi{x;W) {uJi/2{W)cose- kosme)u5^^{x;W) 
+ (1^1/2 (VF) sine + fco cose)u5,e(a;; W) , 
M5^e(x; W) — koUi{x; W) sine + 1/5(2;; W^) cose , 
U5,e{x; W) = koUi{x; W) cose — U5{x; W) sine , 



^This structure of confirms that the deficiency indices of H arc m± = 1. 
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where M5^e(x;W) and U5^e{x;W) are solutions of eq. ^ real-entire in W and 
U5^e{x;W) satisfies boundary conditions ([50| . wc find 



ko 



fl{W) — [/co sine — wi/2(M^) cose] ^ [aji/2(W^) sine + fep cose] 



U5.e{x;W)u5Ay;W), x>y 
U5Ax;W)u5Ay;W), x <y 



the second summand in G{x, y; W) is real for real W — E. 

It is easy to verify that the guiding functional given by (jl7p with U = us^e 
satisfies the properties l)-3) cited in subsec. 13.11 whence it follows that the 
spectra of i?5,e are simple. 

The derivative of the spectral function is given by a {E) — (Trfco)^^ Imil(£'+ 

iO). 

We first consider the case of e = 7r/2 where we have u^^Tr/2{x\ W) — kQUi{x; W) 
and 

<7'{E) ^ (7rfcg)"^Inil](£; + iO) , 
hiW) = giiAay2) + 91 ln(A/fco) - A/2 . 
FoT E ^ p'^ >0, p>0, 2pe^*''/^ we have 

^ ' 2fc2sinh(7r|5i|/2p) " ' 

For E — — < 0, r > 0, A = 2t, and gi > 0, ai/2 = 1 + gi/2T, the function 
il,{E) is finite and real, whence it follows that there are no negative spectrum 
points. 

1/2 



For E — —T^ < 0, r > 0, A = 2t, and gi < 0, ai/2 = 1 — |gi|/2r, we have 



+ 



c 9l ^ _ 2 ^ \g,\ 



{2 + 2nY ^ ko \2 + 2n 

It is easy to see that for the case of e = — 7r/2, we obtain the same results 
for spectrum and eigenfunctions as it must be. 

We thus obtain that for 171 > 0, the spectrum of H^^^^2 is simple, continu- 
ous, and given by spec/fs -1-^/2 ^-nd the set of generalized eigenfunctions 
Ue{x) — \/ a' {E)u^ .^/2{x', E), E >0, form a complete orthonormalized system 
inL2(R+). 

For gi < 0, the spectrum of i?5.±7r/2 is simple and given by specif5_±7r/2 = 
M+ U {£n, n G Z+l and the set of (generalized) eigenfunctions 

UEix) - ^c7'{E)u5,^/2{x;E), E>0, 
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form a complete orthonormalized system in 

We now turn to the case |e| < 7r/2 where we have 

a'{E) = {Trcos^ey^lin[uj5{E + iO)r^ , ujsiW) = fco tan e - Wi/2(W^) . 
For gi < 0, E ^ p"^ > 0, p> 0, X = 2pe-*''/^ we obtain that 

a'iE) = (.cos^e)-Im.-(i.) . , , (51) 

where ujz{E) = A{E) — iB{E). The fmiction B{E) is exphcitly given by 

B{E) ^ I , >0, Vp>0 . (52) 

^ ' 2sinh(7r|gi|/2p) ' ^ ^ 

It follows that for gi < 0, E > 0, the spectrum of H^ ,, is purely continuous. 

For gi > 0, E = p'^ > 0, p > 0, \ = 2pe^^^^^, the derivative of the spectral 
function is also given by eqs. ((5T|) and ((52|) . But in this case, we have B{0) = 
and the limit \imw~^o ^5{W) has to be carefully examined. For small W, we 
have 

UJ5{W) = (tane-taneo)fco- —1^ + 0(1^^) , 

3.91 

taneo = (.91/^0) [ln(.9i/fco) + C - 1] . 

For e eo, the function (j'{E) has a finite limit as i? — >■ 0. But for e = eo and 
small E, we have 



a'(S) = Im + lO)-^ + 0(1) = ^^(5(£;) + O (1) 

TT cos^ eo cos^ eo 



which means that the spectrum of the Hamiltonian H^^ea has an eigenvalue 
E^Q. 

For E — — < 0, T > 0, A = 2t, the function uj^{E) is real. Therefore, 
(j'{E) can differ from zero only at zero points En = En{e) of uj5{E), and cr'{E) 
is represented as 

a'{E) = [-^5(-E„)]"' 5(£^ - En), cusiEn) = 0, io'^{En) < . 

n 

For gi > 0, we have 

uJsiE) = -5iV'(l +gi/2r) - .91 ln(2r/gi) + t + fco(tane - taneo) , 

UJ5{E) ^ ^/\E\- {gi/2)\n\E\+0{l), E ^ -00; a;5(0) = fco (tan e - taneo) . 

For e > eo, the equation ijJ^{E) = has no solution, while for e G (— 7r/2, eo] 
it has a unique solution (e). It is easy to see that 

d,E^~^ (e) = -fco[cJ5 {e[-^^ cos^ e]-i > , 
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so that £;(-) (e) increases monotonically from —oo to when e changes from 
-7r/2 + to eo. 

For < 0, we have 

uj^iE) = I51IVXI/2 - |5i|/2r) + Iffil ln(2r/A;o) + r - g , 
e giC -51 - fcotane . 

Representing the equation ujc^{En) = in the equivalent form 

h{E^) = e, U{E) = |.gi|7/^(l/2- |gi|/2T) + |.gi | ln(2r/fco) + r , 

we can see that: 

a) 

/5(i?) 00, /s ± 0) = ±00 , 

such that in each region of energy (£■„,£„+!), n G (— 1)UZ+, the equation 
i^siEn) = has one solution -E„(e) for any fixed e, |e| < 7r/2, and En{e) 
increases monotonically from £„ + to fn^-i — as e changes from — 7r/2 + 
to 7r/2 — (here, by the definition, £-1 — —00). 

b) For any fixed e, En{e) = —gl/An^ + 0{n^^) as 71 — ^ 00, asymptotically 
coinciding with 



c) The point _E = is an accumulation point of discrete spectrum for gi < 0. 
Note the relation 

lim En-i{e) = lim En{e) = f„, n G Z+ . 

e->7r/2 e-s--7r/2 

The above results can be briefly summarized as follows. 
For gi < 0, the spectrum of H^^e is simple and given by specH^^t = U 
{En < 0, n E (—1) U Z+}. The (generalized) eigenfunctions 



Ue{x) = ^a'{E)u^,,{x-E), E>0 , 

Unix) = [^Lo'^^^iEn]] '^'^ U^,,{x- En), En<0,ne (-1) U Z+ , 

form a complete orthonormalized system in (R+)- 

For gi > 0, the spectrum of iJs^e is simple and given by specH5^e — K+ U 
{E^~'' (e) < 0}. For e G (— 7r/2,eo] the (generalized) eigenfunctions 

Ue{x) = ,Ja'{E)u5,e{x;E), E>0, U{x) = [-uj'^{E'^-^)'\ '^'wg.e (^x;E'^-^^ 

form a complete orthonormalized system in (R+)- For e > eo, the spectrum 
has no negative eigenvalues. 

We note that the above results (for spectrum and eigenfunctions) can be 
extracted from the results in subsec. 13.21 for the case 52 7^ (/i 7^ 1/2). 
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4 Some concluding remarks 



We would like to finish our consideration with a remark about the Kratzer 
potential 1 mentioned in the Introduction. This potential corresponds to a 
particular case of parameters 52 > and gi < 0. It is drown by the thick line in 
the graph of Figure 1. As was already said, the Kratzer potential is extensively 
used to describe the molecular structure and interactions [IF. In such cases, 
the Kratzer potential appears in the radial part of the Schrodinger equation ^ 
and has the form: 



where Df, is the dissociation energy and a is the equilibrium inter-nuclear sep- 
aration. As X goes to zero, V (x) goes to infinity, describing the internuclear 
repulsion and, as x goes to infinity, V (x) goes to zero, describing the decompo- 
sitions of molecules. Putting the potential (|53p in the radial equation 1^ and 
comparing with the Schrodinger equation we have the following identifica- 
tion: 



We can now calculate the value of 52 for real diatomic molecules. Using data 
from [20], even for / = 0, we have 52 = 4.53 x IC* for CO. The parameter 
(72 is of the same order for molecules of NO, O2, I2, and H2. Thus, we can 
see that for the realistic Kratzer potentials, the corresponding radial equations 
have always 52 > 3/4. Thus, the corresponding radial problem belongs to the 
first range described in subsec. 13.11 In this case, there exist only one s.a. radial 
Hamiltonian defined on the natural domain ([5|), functions from this domain have 
asymptotics (fTH]) . 
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